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Highly-symmetric three-band Hubbard Cu-O clusters have peculiar properties when the hole 
number is such that they admit W=0 hole pairs. These are two-hole eigenstates of the on-site 
Hubbard repulsion with eigenvalue 0, get bound by correlation effects when dressed by the interaction 
with the background, and cause superconducting flux quantization. We study the CU5O4 cluster 
by exact diagonalization and show that bound electron pairs of 1 B2 symmetry are obtained at an 
appropriate filling, and quantize flux like the hole pairs. The basic mechanism for pairing in this 
model is the second-order exchange diagram, and an approximate charge conjugation symmetry 
holds between electron and hole pairs. Further, the flux quantization property requires that the 
W=0 pairs of d symmetry have s symmetry couterparts, still with W=0; the former are due to a spin 
fluctuation, while the latter arise from a charge fluctuation mechanism. The simultaneous existence 
of both is an essential property of our model and is required for any model of superconducting d 
pairs. 
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I. INTRODUCTION 

Increasing experimental evidence obtained in several 
cuprate superconductors has convincingly demonstrated 
that the pairing state of these materials has d — wave 
symmetry jlj and that the pairs exist above the critical 
temperature either in the form of superconducting fluc- 
tuations or preformed pairs. The latter aspect is appar- 
ent in the underdoped (normal) region in which a clear 
pseudogap essentially of the same magnitude as the su- 
perconducting gap is measured j|. All these signatures 
put strict constraints to any microscopic model of the 
cuprates. Any theory of the paired state must predict the 
correct symmetry and doping dependence of the binding 
energy of the pair, but the pairing mechanism must be 
doping-indcpendent and robust enough to survive super- 
conducting fluctuations well into the normal state, far 
from optimum doping. 

In BCS theory, the first-order repulsion between like 
charges is overcome by the second order interaction 
with phonons. In high-Tc superconductors the electron- 
phonon interaction is strong and phonons must be ex- 
pected to contribute in an important way to the pair- 
ing interaction, although their task looks harder because 
the repulsion integral U is large (several eV). However, 
the straightforward idea that the hight-T c phenomena 
are just a rescaled version of BCS theory is not granted. 
The role of phonons may be important, but is differ- 
ent, and some other ingredient is essential. First, due 
to the planar d v symmetry of these materials there is 
actually no repulsion barrier to overcome. In a series of 
papers ||-[| we have introduced the two- hole X B2 eigen- 
states with zero Coulomb on-site repulsion (the so called 
W=0 eigenstates). This is based on a standard two- 
dimensional Hubbard description of the particle corre- 
lations in the copper-oxide plane. The model has been 
numerically implemented in real space on highly symmet- 
ric Cu-0 clusters containing up to 21 atoms and there the 
W=0 pairs arise from degenerate (x,y) levels. We have 
shown that when the first order interaction vanishes the 
net effect of the interaction is attractive in the physi- 
cal parameter range. We used the smallest hole number 
(n=4 holes) necessary to form the W=0 pair on the de- 
generate e(x, y) one-hole level of the clusters. Also, we 
have shown that although U is not small the dressing of 
W=0 pairs is perturbative; by exact diagonalization and 
diagrammatic analysis we demonstrated that the model 
leads to pair attraction of the order of tens of meV aris- 
ing from virtual spin-flip excitations. Thus, rather than 
proposing a model of pairing, we have pointed out that 
the most standard description of electron correlations in 
the Cu-0 plane already leads to pairing in clusters when 
the symmmetry of the problem is fully allowed. Remark- 
ably, the off-site repulsive interactions, when included, 
tend to enhance the effect somewhat Q, so we devote 
the present study to the on-site interaction effects for 
simplicity. The mechanism we are considering is only 



a part of the story, but it seems to be the most pecu- 
liar part, being related to nothing but the C4 U symme- 
try. As such, the mechanism itself is doping-indepcndent 
(although its effects do depend on doping). For similar 
reasons here we wish to make abstraction from phonon ef- 
fects to see how far the simplified mechanism can account 
for reality by itself. We believe that a mechanism which 
predictably gets attraction out of repulsion is by itself of 
theoretical interest. Superconductivity is a much more 
complex problem than pairing, and the thermodynamic 
limit of the present model is currently under considera- 
tion; however, we pointed out ]5| that the magnetic flux 
is quantized by our pairs in the same way as it is in the 
type II superconductors. 

In the present paper we extend the analysis of Ref. (|] 
by diagonalizing the CU5O4 cluster with increasing num- 
ber n of holes. Eventually, we find electron pairing, again 
with a binding energy A(n) of a few tenth of eV in the 
physical parameter space. Electron pairing is actually 
realized (6), e.g., in the T' structure of (Nd,Ce)2Cu04 
exhibiting superconductivity. The T' structure of this 
compound is different from the T structure of La2Cu04, 
but is still characterized by CuC>2 planes 0.Our main 
point here is that electron pairs and hole pairs are re- 
lated by an approximate charge conjugation symmetry 
and the very same basic mechanism or diagram is oper- 
ating in both cases. Further, we demonstrate how the 
two different symmetries (Ai and B2) of W=0 singlet 
pairs allowed by our theory produce the superconduct- 
ing flux quantization phenomenon, both for electron and 
hole pairs, in the respective doping regimes, and are thus 
both necessary for a serious proposal of a mechanism of 
pairing. 

II. CHARGE CONJUGATION 

By a canonical transformation from holes to electrons, 
the three-band Hubbard Hamiltonian we considered pre- 
viously H becomes: 

H = ^(2e, + Ui) - Y^i £ i + U z )a+a la 

i icr 

- ^ hjal^aja +22, u i n i+ n i--> (!) 

<i,j>cr i 

with a^ T =Ci cr and ai (T =c^ T creating or destructing an elec- 
tron of spin <j= f , J, at site i, respectively. Thus, in the 
electron representation the signs of the site and hop- 
ping integrals are reversed. The single-electron energy 
levels of the CU5O4 cluster have been computed in the 
physical range of parameters (in eV) used before, i.e., 
Ud=5.3,U p = 6, t=1.3,£ £ 2=0,£ p =3.5 and the results are 
displayed in Fig.l as a function of the oxygen-oxygen 
hopping integral t ox . With respect to the hole case [||, 
the sequence of levels is such that bonding and antibond- 
ing states are interchanged. According to symmetry ar- 
guments j| , singlet pair eigenstates of the Hamiltonian 
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H with zero eigenvalue of the interaction part and : B2 
symmetry exist when two particles sit on a degenerate 
e(x, y) level. In the hole representation the antibonding 
state of e symmetry contains two holes when 14 holes 
are present in the cluster. Its charge-conjugate counter- 
part (since the dimensionality of the one-body basis is 
18) corresponds to a total of 4 electrons in the cluster, 
two of them filling the lowest a,\ bonding state and the 
other two sitting on the next e(x,y) level of Fig.l. 
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high-T c cuprate in which charged carriers are electrons. 
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FIG. 1. One-electron energy levels of the C115O4 cluster versus 
ton. Parameters are given in the text. The levels are labelled 
according to the representations of the C4,, Group. 



In the symmetric CU5O4 cluster we previously found 
H that the exact diagonalization of the Hamiltonian ma- 
trix with four holes in a wide range of negative t ox val- 
ues for the 1 B2 singlet gives a negative pairing energy 
A(4), where A(n) =E(n)+E(n-2)-2 E(n-l), and E(n) is 
the ground state energy with n holes. Therefore one ex- 
pects that A(14) is negative while A(12), which corre- 
sponds to the complete filling of the upper e(x,y) level 
with 1 Aisymmetry, is highly repulsive. Typical energies 
A computed with t ox =- 0.2 eV and different hole num- 
bers are listed in Table 1. 

The t ox dependence of A computed using an enhanced 
Lanczos diagonalization routine is shown in Fig. 2 and 
Fig.3 for A(4) and A(14) and for A(10) and A(12), re- 
spectively. Filling the levels with n holes, pairing occurs 
when the uppermost pair is a W=0 singlet. It is appar- 
ent that a considerable degree of symmetry under e-h 
exchange exists in the cluster and that the pairing in- 
teraction is not restricted to a single value of doping. 
The electron- hole symmetry is not exact, because the 
one-electron energy level spectrum of Fig. 1 is not in- 
variant when up is exchanged with down and left with 
right (t ox — > —t ox and E — > —E); however one could 
have predicted the negative A and its magnitude by the 
approximate symmetry. In real systems electron super- 
conductivity has been first reported Q in the Nd2CuC>4 
doped with Ce with a maximum T c =24 K for optimal 
doping concentration x=0.15, similarly to more common 
hole superconductors. Nd2- :r Ce :r Cu04 is an example of 
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FIG. 2. Pairing energy A(n)=E(n)+E(n-2)-2 E(n-l) of the 
CU5O4 cluster with n=4 and n=14 holes. 




FIG. 3. A(n) as a function of t ox for n=10 and n=12. 



The spin-flip diagram responsible for the pairing has 
been identified before ||. The shift of the pair energy 
was found, after considerable algebra, in terms of the 
one-body site amplitudes. When A (4) is computed, the 
background holes are just 2, and occupy an orbital of a± 
symmetry, which we denote by a. Empty states (anti- 
bonding and nonbonding) of the same symmetry will be 
denoted by a', and we shall write b for the empty states 
of bi symmetry. To second-order one obtains for A the 
following expression: 



A = -2 



W(a, b, x, xf 



. b 



e 6 



E 



W(a, a' , x, x) 



where W{k, l,m, n) 



: (k + l-\Y J U l n i 



(2) 



Mi— |m_i_n_) and 



the sums run only over the one-body states of a and 
b symmetry since the contribution of the e(x,y) orbitals 
vanishes. The final sign of A is determined by the rela- 
tive weight of the virtual excitations to the empty states 
of different point symmetry. 
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It is important to check the range of validity of the 
perturbative expansion by comparing the results of Eq.2 
with the exact diagonalization values of A. We construct 
the normalized hole states of symmetry a, b, e(x, y) whose 
site amplitude and energy are listed in Table II. Making 
use of the above standard values of the parameters we 
calculate the accuracy 8 — 2 ^ ej; , ^ appr for n=4 as a 
function of negative t Q 



A ex +A a 

's, i.e. in the region of pairing. 
The accuracy is good (about 3%) for small U values ( 
< 0.05 t ) and raises steeply to approximately 30% for 
larger values of U/t with a smooth dependence on t ox . 
However, for t ox =0 where the pairing is strongest for the 
1 B 2 singlet, the accuracy is around 7% up to U/t «1. 

The values Ud = 5.3 eV, U p = 6. eV differ appreciably 
from other literature estimates ||, and must depend on 
the compound and doping. For I^CuCU, J7 P =4 eV and 
Ud=l0.5 eV have been recommended [||. None of the 
above results depends qualitatively on the precise value of 
the model parameters, since ours is basically a symmetry 
argument. 



III. FLUX QUANTIZATION AND PAIR 
SYMMETRY 

Consider the ground-state energy E(0) of a two- di- 
mensional system as a function of the magnetic flux 
through it. This is definitely a most compelling way of 
testing its superconductivity, since type II superconduc- 
tors quantize the flux for integer and half-integer multi- 
ples of = — . If one inserts a flux tube in a loop formed 
by closing the path on the external Cu of the CU5O4 clus- 
ter, one finds || that the calculated ground state energy 
E(4) has a second minimum at 0o/2 well separated from 
a barrier from the 0=0 minimum. The point symmetry 
of the wavefunction changes from 1 B2(x 2 — y 2 ) at 0=0 to 
1 Ai(s 2 +y 2 ) at 0=0n/2. Owing to the symmetry under- 
charge conjugation and to the attractive interaction of 
the E(14) state, one expects to find the specular situa- 
tion in the electron case. That this is indeed the case 
it is apparent from Fig. 4 where the response functions 
R=[E(0)-E(O)]/tc; of the hole- and electron-doped clus- 
ter are compared. Here t^ is an (in principle, infinites- 
imal) hopping energy between the external Cu's, which 
is needed to close the loop around the tube. Due to the 
small size of the system the barrier preventing the con- 
densation of pairs is finite but of comparable height in 
both cases. The symmetry of the states is the same as in 
the hole case, i.e., 1 B2 for 0=0 and x Ai for 0=0n/2. 
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FIG. 4. Response function R versus the normalized flux for 
the CU5O4 cluster with n=4 and n=14 holes. 



However, the mechanism of interaction between the 
holes (electrons) is basically different for 1 B2 and 1 A\ 
W=0 pairs. Consider again the interesting case t oa: =0 
for which all the oxygens are equivalent and the point 
symmetry Group is S4. Within the space of B2 sym- 
metry spanned by the |xy> and |yx> two-body states 
the contribution to the scattering amplitude is zero to 
first-order, but finite and negative to second order. Thus 
the anomalous propagator A s f ps (t^-L^ sca tters the 
|y + x_ > state out of the |x + y_ > state where A^ 2 ) is 
given in Eq.2 and e is the energy of the unperturbed 
e(x, y) states. The interaction is attractive for the low- 
est state of 1 B2 symmetry and involves a second order 
spin-flip fluctuation. 

Consider now the basis of the degenerate |xx>+| yy> 
states and the |bb> states all having : Ai symmetry. The 
eigenvalue problem requires the diagonalization of the 
2x2 matrix: 



< bb\W\bb > 
< bb\W\x 2 



< bb\W\x 2 + y 2 > 
y 2 > < x 2 + y 2 \W\x 2 + y 2 > 



Up Up_ 

U P Up_ 
2V2 2 



(3) 



where U p is the on-site repulsion at the oxygen site. 
The lowest eigenvalue is and the W=0 pair is ip — 

The upper eigenvalue is 3U/4 



§!&&)• 



II- 2 - 



y 2 > 



and the eigenfunction ip = y + y ||a; 2 +y 2 > . Thus 
we have a first-order mechanism involving a charge fluc- 
tuation. 

By looking at Fig. 4 we realize that with increasing the 
flux 0, the X B 2 component of the ground state is quickly 
destabilized by the vector potential, while the : Ai com- 
ponent decreases in energy and eventually becomes ener- 
getically favored in presence of the trapped flux. A simi- 
lar change of symmetry occurs in conventional supercon- 
ductors where the Cooper wavefunction has s symmetry 
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fL(| . The Hubbard large U system is stabilized by the in- 
trinsic charge fluctuation. Since the mechanism is robust 
(first order) it is conceivable that these low lying pairing 
states can survive in the vortex core of the underdoped 
phase |Q. The 1 Ai symmetry of the ground state in the 
presence of the half flux quantum is required by general 
symmetry principles. Indeed, the vector potential lowers 
the symmetry from C^ v to its abelian subgroup Z4 , which 
contains only the rotations. In Z^Bi and 1 B2inerge into 
the same irreducibe representation and so the vector po- 
tential mixes them, giving rise to pairs. This does 
not apply to the 1 Ai state, which remains single and 
W=0. With increasing the flux, the component of 
the ground state increases its energy, until it crosses the 
1 Ai pair energy. The minimum at half flux quantum cor- 
responds again to a A <0 situation. It is gratifying that 
the present model allows for W=0 pairs of both symme- 
tries, because this is the only way the superconducting 
flux quantization works. 

IV. CONCLUSIONS 

We have investigated the symmetry properties of the 
W=0 pairs under charge conjugation by performing nu- 
merical diagonalizations of the CU5O4 cluster doped with 
a variable number of holes n ranging from to 14. We 
find attractive interaction of comparable strength when- 
ever the bonding and antibonding degenerate levels of 
x — y symmetry are half-filled. All the other fillings lead 
to strong repulsion. The 14-hole case corresponds to the 
doping with 4 electrons and thus to electron supercon- 
ductivity. In all cases the attraction is provided by the 
spin-flip fluctuation between |x+y_ > and |y+x_ > pairs. 
The n=4 and 14 hole 1 B2 states quantize the supercon- 
ducting flux and the ground states in presence of zero 
and one fluxon have different symmetry as in BCS su- 
perconductors. 
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TABLE CAPTIONS 

TABLE I. A values of the CU5O4 cluster vs the number 
of holes n computed with the standard set of parameters 
given in the text and t ox =-0.2 eV. 

TABLE II. Symmetry, site amplitude and energy of 
one-hole levels for the CU5O4 cluster. In the first row we 
report the Cartesian coordinate of each site in units of the 
Cu-0 distance. The first column shows the orbitals; clq 
is the nonbonding level belonging to the Irrep A\ . Minus 
and plus superscripts refer to bonding and antibonding 
levels, respectively. 



TABLE I 



n holes 


A(mcV) 


4 


-15.9 


10 


141.2 


12 


914.7 


14 


-22.4 
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TABLE II 
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